Harmonium in intense laser fields: excitation, 
absorption, and transparency 

O. Kidun and D. Bauer 

Max-Planck-Institut fur Kernphysik, Postfach 10 39 80, 69029 Heidelberg, Germany 

Abstract. It is known that the dynamics of two (Coulomb-interacting) nonrelativis- 
tic electrons confined by a parabolic potential and driven by a classical, intense laser 
field (in dipole approximation) is exactly soluble. We calculate the time-dependent 
population of the harmonic oscillator states and the energy absorbed from the laser. It 
turns out that the key entity on which all observables sensitively depends is the mod- 
ulus square of the Fourier-transformed vector potential of the laser field, evaluated at 
the harmonic oscillator frequency. The system is transparent to laser field configura- 
tions for which this entity vanishes. We discuss the Poisson statistics behavior of the 
transition probabilities and analyze the conditions for the complete survival and full 
depletion of the initial state. 



1. Introduction 

The parabolic well is a model of great practical importance, as it approximates 
any arbitrary potential close to equilibrium. In nanotechnology, potentials of simple 
shape such as quantum dots are often well approximated by parabolic potentials. 
Superpositions of harmonic oscillators are also used to describe continua, environments, 
and field modes. 

Correlations, for instance, those introduced by the Coulomb interaction between 
electrons, attract growing attention. Correlated electron dynamics can nowadays be 
analyzed experimentally with great precision using "reaction microscopes" [I] . Two 
electrons confined in a parabolic potential, sometimes called "harmonium" or the 
"Hooke's atom" version of helium [2], is of great interest since the problem represents 
one of the few exactly soluble problems involving correlation [31 H] . As such harmonium 
serves as a testing ground for, e.g. new exchange- correlation potentials in density 
functional theory. 

In the field of strong laser atom interaction, exactly soluble models are rare as 
well, in particular as it comes to electron correlation. Fortunately, the two-electron 
wavefunction of the laser-driven harmonium can still be expressed in analytical form as 
long as the laser field is treated classically and in dipole approximation [5]. 

The outline of the paper is as follows. In Sec. [2] we derive the exact solution of the 
time-dependent Schrodinger equation of two interacting electrons confined in a parabolic 
potential and driven by a classical (laser) field (in dipole approximation). In Sec. [3] the 
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probabilities for two-electron excitations from an arbitrary state m to a state n are 
calculated. The latter are the basis for the discussion of optimal energy absorption and 
transparency in Sec. HI Finally, we conclude in Sec. [5j 



2. Exact states of two-electrons in a parabolic well driven by a laser field 

In the absence of the laser field the Hamiltonian of harmonium read 

H{IU r2) = _£v. - + + j^i + ^ (i) 

with it) the harmonic oscillator frequency and m the mass of the electron. The 
Hamiltonian (TjQ) is separable when written in the center-of-mass (CM) and relative 
coordinates R = |(ri + r 2 ) and r = r x — r 2 , respectively: 

H = H R + H r = f-^V| + ™ 2 R 2 ) + {--Vl + ^Vr 2 + ^ . (2) 
y 4m ) \ m 4 \r\J 

As a consequence, the eigenfunctions $(ri, r 2 ) are products of the eigenfunctions of the 
CM and relative Hamiltonians Hr and H r , i.e. with if^ \<f>) = Er \(f>) and H r |£) = i? r |£) 
follows 

$( ri , r 2 )e-* E * = $(R, r)e-^ E « +Er)t = 0(R) e -^«* • ^(r) e -^ rt . (3) 

Due to the spherical symmetry of both Hr and H r the eigenfunctions are conveniently 
expressed in the form (j) NLM (R) = Unl r {K) Y lm (Q,x) and ^ m (r) = ^j^-Y em (6, tp). 

The CM radial wavefunction and energy are readily obtained from the standard 
problem of the 3D harmonic oscillator to give Unl{R) = CNLR L+1 H N (R)e~ m,JjR2 ^ h 
and E R = hto(2N + L + §). Hn^F^-N, L + l,muR 2 /h) is the iV-th order Hermite 
polynomial, Cnl is the normalization constant. The radial wavefunction u n e(r) has the 
closed analytical form (unnormalized) u n e(r)=r e+1 e~ mujr2 ^ 4h J2k=o A k r k only for certain 
oscillator frequencies [3] • For other frequencies, it has to be evaluated numerically. Once 
it is found, the eigenenergy has the simple form E r = Tiuj{n + t + |). 

Adding the interaction with an electromagnetic field affects only the CM 
Hamiltonian so that the problem still separates. Having in mind the interaction with 
near infrared or visible laser pulses, we adopt the dipole approximation and obtain 



h 2 _r, m r. e 2 \ ( 1 



H= \ —Vt+'—u 2 r 2 + 




2e 

-iW R +— A(t) 

c 



i 2 



+ muj 2 TL 2 ) (4) 



m ' r ' 4 |r| J ' \2 ■ 2m 

with A(t) the vector potential and c the light velocity. The electric field of the laser is 
given by E = — d t A while the magnetic field is neglected^ in dipole approximation. The 
total wave function reads 

*(R,r,t) = V(R,*)-e(r)e-^ (5) 

with £(r) the same as before. 



| Spin is omitted for brevity. 

§ Two parabolically confined electrons in constant magnetic fields were studied in [6] . 
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The solution of the time-dependent Schrodinger equation governing the CM motion 
with twice the electron mass \i = 2m and charge e = 2e 



-ihV R + -A(t) 



+ -^ 2 R 2 U(R,t) 







(6) 



is known [7J [8]. In the case of a linearly polarized laser field A(t) = A(t)e z one has 

^(R, t) = U(X)e-^ Ext U(Y)e-* EYt U(Z, t), (7) 
with the two unaffected harmonic oscillator eigenfunctions U(X) = 



-iiwX 1 /2ft 



U(Y) = C Ny H Ny (X)q ^ y l 2n and the respective eigenen- 



ergies Ex = hcu(Nx + §), Ey=huj{Ny + §)■ Instead, the eigenfunction U^ Z (Z) becomes 
dressed by the laser field and reads 
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(8) 



The laser field 



The normalization constants are given by Ck = 
is turned on at time t = t . Before, the system is assumed to be in an eigenstate 
determined by the quantum numbers N x , Ny, and N Z - The functions f(t) and g(t) are 
given by 



f(t) = cos tot / A(t) sinturdr — sin cut / A(r) cosc^rdr , 

Jto Jto 

g(t) = sin cut / A(r) sincurdr + cos cut / A(r) coscurdr . 

J to Jto 



(9) 
(10) 



^From the structure of (jSJ) one can infer (i) that the center of the CM wave 
packet describes motion along a trajectory R(£) = [0, 0, — g(t)] where, in fact, g(t) 
is proportional to the excursion of a driven, classical harmonic oscillator, and (ii) that 
the time-dependent solution U (Z,t) represents a so-called coherent wave packet, i.e. a 
state of minimum uncertainty Ap AZ = h/2, equally distributed over the momentum 
p and the spatial coordinate Z (see, e.g. [9]). 



3. Photoexcitation of the electron pair 

Let us consider the photoexcitation of an electron pair confined in a parabolic well 
(e.g. two electrons occupying the low- lying states of a quantum dot). The transition 
amplitude is given by the overlap of the exact time-dependent two-electron wave function 
ty(ri,T2,t) with the asymptotic stationary solution $(ri,r 2 ) when the field is switched 
off: 

T(t) = <$|tf(t)> = (cPnlm I ll>N.N,Njit)) = (U Kz | U Nz {t)) . (11) 
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Here we made use of the fact that neither the sub-problem of relative motion of the two 
electrons nor the X and Y components of the CM motion are affected by the laser field. 
What is left is the probability to find the CM quasiparticle occupying a stationary state 
U N fin(Z) of the free oscillator after the action of the laser. 

Let us first discuss the case where we start from the ground state, i.e. N x = N Y = 
Ng 1 = 0, (R\il;(t)) = i[) 00 o(Il,t). The squared modulus of the corresponding transition 
amplitude to some final state with the quantum number N z n = n reads 

r2,_, „2,_, A ' 2 (r) 



r (o-*»)( t ) 



u u n 



exp < mujt + i 



/oo 
dZH n 
oo 




Z\ exr 




(12) 



One can omit all purely time-dependent imaginary exponents, producing unity, and 
reduce the integral to the table form [10J by change of variables 

2 



ir(t)r 



2fihc 



/oo 
dx H n (x)e~ 
oo 



(x-vY 



(13) 



with x 



^Z and y 



As a result, the probabilities to find the CM of the electron pair in the ground (that 
is the survival probability) or in the n-th excited state are 

(14) 

2 O n 



P^°\t) 
p(0-n)( f ) 

a(t) 



T^°\t) 



nl 



e 2 uj 



2^hc 2 



9 2 {t) + f{t) 



(15) 
(16) 



where we introduce the important parameter a(t), which depends on the laser field and 
the oscillator frequency [cf. eqs. © and tTTO]) ] . 

The formula for the photoexcitation probability (TT5l) describes the distribution of 
the electronic wave packet among the oscillator states as a function of time and the EM 
field parameter. It has the form of a Poisson distribution, which is characteristic for 
stochastic processes where the system absorbs the first amount of energy from the field 
independently from the possibility of the absorption of the second portion and so on. 

In Fig. [1] we illustrate the above expression for p(°^ n ) for a laser pulse of frequency 
Mi = 2 atomic units (au) with a Gaussian envelope and halfwidth t p — 1 au, A(t) = 
sin uit exp (-|(t/t p ) 2 ) (dashed curve). The survival probability P °(t), changing from 
unity to some final value, has one maximum during the pulse due to the oscillation of 
the CM wave packet inside the potential and the associated maximum overlap with the 
ground state as the center of the wave packet passes through Z = 0. The transition 
probabilities of the excited states behave in the opposite way, i.e. they are minimum at 
these times. 

In the same way as for the ground state excitations [cf. eqs. (1111 — [T5]l ]. the general 
expression for the transition probability from an arbitrary initial m th oscillator state 
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(t < t ) to the final n th state can be found; it reads 



p(— )(t) = 2 F Q 2 (-m,-n;—) — e" CT , (17) 



a J n\m\ 

with the hypergeometric function 

- -i . min(n,m) . . -, 

_, / 1 \ ±A ni ml 1 

2 F -m, -n; — = }^ -rr t-t - . 

V a J [n — k)l (m — k)l ki \ a 

Equation (fTTj) circumscribes the time dependent population of the two-electron 
excited states following the laser pulse. As the Poisson distribution ffl~5]) characterizes the 
stochastic absorption of n photons, the distribution ( 1171) is responsible for two processes, 
namely the independent emission of m — k photons and the (re) absorption of n — k 
photons. The polynomial sum (fl8l) in the expression accounts for the transitions through 
all the allowed intermediate states k. 

Of particular interest are the asymptotic excitation probabilities lim^ +00 p( m ^™) (t) 
since they are experimentally accessible observables. In Fig. [2] the distribution p( m_ * n ) 
is portrayed as a 2D intensity plot vs the (analytically continued) excitation quantum 
number n and the field-strength parameter a. The latter is proportional to the laser 
intensity (see the discussion in the next section). Plots (a) and (b) correspond to the 
different initial states with m = (ground state) and m = 3 accordingly. For m = the 
quantum number of the maximally populated final state n* is governed by the equation 

d p(O-n) = d / 1 = Q) (lg) 



dn dn \n 

which leads to the simple relation 



a = e A( "* +1) = n* + - + o (lCr 1 ^"*) . (20) 



A( n ) = Tn^i x ) i s ^ e Euler polygamma function. 



dx " ' x=0 



4. Energy absorption 

Let us now calculate the energy absorbed from the laser field by the two-electron 
system initially prepared in the ground state E ahs = E{t-^ oo) — E(t = t ), where 
E{t) = (ili (t)\H(t)\ifa){t)) with H(t) the total Hamiltoniaijjjl 



2, ,2 



e u 



2/ic 2 



A(t) A 2 (t)\ fiu 



(Mt)\H(t)\Mt)} = At) + f(t) +2f(t)^ + — + —, (21) 



, ^ 2 

E a bs = hm ^ 



uo uo 2 } 2 



g 2 {t) + f{t) = huo lim a(t). (22) 



This expression for the absorbed energy is obtained independently of the initial state, 
i.e. using in eq. ( [211) the wave function \ip m {t)) (starting from an initial state m) also 
leads to the result fl22|) . The absorbed energy depends on lim 4 ^ 00 cr(t) only but not on 

|| We disregard the possible spontaneous decay of the excited system. 
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the initial state or the population over the final states of the system. The parameter of 
central importance thus clearly is a(t — > oo) = o~ Q , which can be recast into 



2, , A 



o n = lim 



e oj 



t— >oo 



e u) 



g\t) + f\t) =-— \i m \g(t)+if(t)\ 



t— >oo 



6 2 U 



2 fihc 2 



+oo 



dt'A{t'] 



e 2 cu 



\A{u>)\', (23) 



2/iftc 2 

showing that the spectral power of the incoming laser field evaluated at the harmonic 
oscillator frequency determines the energy absorption E a ^ s . 

The expression for the absorbed energy ( 1221) can be alternatively derived via the 
average number (n) of absorbed photons calculated using the photon statistics: 

oo n 

E ahs = Kuj in) = %uj V n e~ a ° = huoa . (24) 
The essential fact is that the field strength parameter contains the cumulative mass 

2 

and the charge of the CM subsystem with different powers cr ~ — . Therefore, the 
observables in the single and in the double-electron excitation differ from each other. In 
particular, the energy absorption from the same laser pulse by a (2e)-system is twice as 
large as in the one electron case. 

In the generalized many-particle problem the parabolic potential contains an 
arbitrary number of electrons N. As in the (2e)-case the iV-electron CM subhamiltonian 
accounting for the EM field can be separated from the field-independent relative motion 
of the electrons. Then, all the formulae describing two-electron transitions are valid for 
the many-electron case up to the effective charge and mass. Hence the energy absorption 
changes linearly with the number of electrons in the well, E^' /E^J = N. Given this 
property, the number of active electrons in a harmonic system may be experimentally 
accessible. 

Let us now discuss the two extreme cases of complete survival and full depletion 
of the initial state, that is its maximum and minimum survival probability. If the 
survival probability at the detection time equals unity, p( rn ^ m )(t — > oo) — > 1, the 
quantum system does not absorb energy from the laser field and becomes transparent 
for the corresponding pulse. This is strictly realized at the zero value of the field- 
strength parameter a a = in the exponent in ffT?|) . Hence, according (1231 . the absence 
of the oscillator frequency in the spectrum of the vector potential, A(uj) = 0, precludes 
absorption, [cf. eq. (122]) ]. Note, that this statement is correct for any field intensity and 
not just in a first order perturbative treatment as for nonparabolic confinements. 

As an example, consider a laser field consisting of two consecutive Gaussians shifted 
in time by the interval 2a: A(t) = sxa[ui(t + a)]e _ 5[(*+ a )/*pl + sm[ui(t — a)]e~5t(*~ a )/*p] . 
Setting e.g., the laser frequency coi = 1 and the halfwidth of pulses t p = y/2/2, we 
aim at finding a relation between the time delay 2a and the oscillator frequency u, for 



which = A(uj; a) = —i^n/y/e cosuja sinh f e~ w / 4 . The result is ua = ir(k + 1/2) 
with k an integer. This is, in fact, a condition for the destructive interference (in time) 
of the two waves within the characteristic time-domain 27r/u of the system. In terms 
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of the excitation picture the explanation is the following: the electronic wave packet, 
distributed over the excited levels by the first pulse, may be assembled back to the initial 
state by the second pulse, if the phase difference between the pulses is chosen properly. 

The opposite case amounts to the minimum survival probability, i.e. the maximum 
depletion of the initial state. On account of e~°"° — ► this is the case if the spectral 
power at the transition energy (coinciding with the oscillator frequency uj) is infinitely 
large. It is equivalent to the presence of an infinite, monochromatic laser pulse 
A(t) = A Q e~ tuJlt with A{uS) ~ A 5(uj — uji). Such a pulse resonantly depletes, for 
instance, the ground state. For a finite pulse, e.g., A(t) = smu[t e~2(*/'p) 2 with a power 
spectrum E^ s ~ \A(uj)\ 2 = 27itpSmh 2 (uJUitp)e~^ 2+UJ ^ t p ^ 0, there is always nonvanishing 
absorption. However, only for the resonance uj\ = uj the absorption exists in the limit 
of large pulse durations, t p — > oo. 

In addition, if the initial state differs from the ground state, its full depletion 
occurs at such values of a Q satisfying a™ - 2-^o( — m , ~ n \ ~ jr) — 0. For example, the 
well prepared state with m = 1 will be exhausted and distributed over other levels 
by the laser field of the strength o Q = 1 [cf. eq. ffT7|) ]: P^ 1 ) = (er — l) 2 e~°°. Similar 
relations can be obtained for other higher order initial states, e.g. o a = 2± \J~2 for m — 2, 

a Q - 3 - a - | = 2a Q - 6 + a - f ± iv^ (a - |) = 0, a = (3 + 3^) for m = 3 and 
so on. 

5. Conclusion 

In this paper we considered the exactly solvable problem of the excitation of two 
interacting electrons, confined in a parabolic well in the presence of a laser field. The 
expression for the time- dependent population of the excited two-electron states was 
derived. The physical features of the process, such as the Poisson statistics of the 
transition probabilities and the linear dependence of the energy absorption on the laser 
power, were analyzed, in particular, the conditions for the complete survival and full 
depletion of the initial state. Our approach can be used for the testing of approximate 
theories dealing with few-body multiphoton excitations. 
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Figure 1. The time evolution of the survival probability in the ground state, n = 
(black curve), and the excitation probabilities to the states with the quantum numbers 
n= 1, 2, 3, 4 (green, blue, orange and olive curves, respectively) is shown together with 
the laser pulse (red dashed line) A(t) = sin. u> it e - '*/*") / 2 . The laser frequency, the 
half- width of the pulse and the oscillator frequency were chosen uji = 2 au, t p = 1 au, 
and ui = 1 au respectively. For better visibility, the fast decrease of the probabilities 
with the quantum number p(°^ n ) ~ i, jg compensated by multiplication with the 
factors given in brackets. 
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Figure 2. The asymptotic distribution p( m ^ n )(t — > oo) of the two-electron wave 
packet over the excited states, i.e. their population vs. the (analytically continued) 
quantum number and the field intensity is illustrated as a 2D intensity plot. Figures 
(a) and (b) correspond to the different initial states with m — (ground state) and 
rn = 3, respectively. The color scale of the population probability is common for both 
figures. The red line in (a) [determined by eq. ([20]) ] denotes the position (the quantum 
number) of the maximally populated state for the given field-intensity parameter a D 
(c.f. eq.$2M). 



